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Abstract 

Modern classical computing devices, except of simplest calculators, have von 
Neumann architecture, i.e., a part of the memory is used for the program and 
a part for the data. It is likely, that anàlogues of such architecture are also de- 
sirable for the future appUcations in quantum computing, Communications and 
control. It is also interesting for the modern theoretical research in the quan- 
tum Information science and raises challenging qüestions about an experimental 
assessment of such a programmable models. Together with some progress in 
the given direction, such ideas encounter specific problems arising from the very 
essence of quantum laws. Currently are known two different ways to overcome 
such problems, sometime denoted as a stochastic and deterministic approach. 
The presented paper is devoted to the second one, that is also may be called the 
programmable quantum networks with pure states. 

In the paper are discussed bàsic principies and theoretical models that can be 
used for the design of such nano-devices, e.g., the conditional quantum dynam- 
ics, the Nielsen-Chuang "no-programming theorem," the idea of deterministic 
and stochastic quantum gates arrays. Both programmable quantum networks 
with finite registers and hybrid models with continuous quantum variables are 
considered. As a bàsic model for the universal programmable quantum net- 
work with pure states and finite program register is chosen a "Control-Shift" 
quantum processor architecture with three buses introduced in earlier works. It 
is shown also, that quantum cellular automata approach to the construction of 
an universal programmable quantum computer often may be considered as the 
particular case of such design. 

1 Introduction 

There are two almost independent ways of the classification of the quantum computa- 
tional networks. In this paper and in many other works about the theory of quantum 
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computations are used models with quantum networks acting on pure states repre- 
sented as vectors in Hilbert spaces [1]. It is also possible to consider mixed states 
and to use density matrices [2]. Such methods are especially useful in theory of open 
quantum systems, but for networks described in this paper it is enough to consider 
only pure states. 

On the other hand, the structure of a quantum network may be specified using few 
different levels, which also have some analogue with the classical case. The simplest case 
— is a fixed network for the resolution of a particular task, e.g., the quantum network 
for Shor's factoring algorithm [3]. In the quantum case there are three different steps: 
the preparation of initial state, the quantum evolution described by given quantum 
network, and the measurement (Fig. la). 

More difíicult level — is a network with tuneable and interchangeable elements for 
managing with different tasks (Fig. Ib). It is common also for most experiments. Here 
is also present the initialization and the read-out (measurement), but between them 
may be considered different sequences of operations. Already in the earliest works 
[4] was raised a question about an universal set of such elementary operations, the 
quantum gates. 
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Figure 1: Usual schemes of quantum networks. a) The simplest quantum network with 
a fixed structure. b) The network with a set of gates for external control. 

Even the second kind of network still rather resembles a quantum "calculator" 
Fig. Ib than a computer, because it is controlled via external manipulations instead 
of doing some program. It is especially essential for the quantum case, because such a 
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control here is described as a classical process and so there is an additional difficulty 
for the unified description of the whole system. 
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Figure 2: The programmable quantum network. 

It is possible to consider a next level — the programmable quantum network with 
two subsystems: the data and the program Fig. 2. It is similar with von Neumann 
architecture of the classical computers [5]. 

With the usual notation of the quantum computation the programmable quantum 
networks U (with pure states) may be described as: 

(j{\D)\n)) = {un\D))\n). (1) 

Here \D) and |i7) are states of data and program registers before operation. After 
the application of the fixed unitary operator U, the state of the data register may 
be described as Un\D), i.e., some operator Un is applied to the data state and it 
depends on the state |i7) of the program register. State of the program register after 
the operation U is not changed (a more general case is represented in Eq. (1') below 
on the page 5). 

It was found [6] that for the pure states Eq. (1) may be vàlid only if all different 
states \n) of the program register are orthogonal. E.g., if there is some program (state) 
|i7i) for the implementation of an operator Ui and we need to implement another 
operator U2 using some program |i72), then: 



(Í7i|i72) = 0. 



(2) 
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The Eq. (2) was dcrivcd in [6] from Eq. (1) and unitarity of U (see Sec. 2.2). 

It was considcrcd in [6] as somc kind of the no-go theorem for universality of a 
programmable quantum network, because in such a case dimension of the Hilbert 
space for the program register must be equal to the number of diíferent programs, i.e., 
unitary operators, but for the exactly universal quantum computer there are infinite 
number of such programs. 

There are a few important no-go results in the theory of quantum computing and 
most known — is a no-cloning theorem [7]. The no-cloning theorem may be deduced 
from linearity of quantum mechanics, and the just mentioned "no-programming" prob- 
lem even more subtle, e.g., exists a linear non-unitary operator satisfying Eq. (1) with 
dimension of the program register only in two times bigger than the data onc [8] (see 
Sec. 2.3). The problem with infinite dimension of the program register in the univer- 
sal programmable quantum network appears due to the orthogonality of the diíferent 
program states Eq. (2) derived from unitarity of U. 

The no-programming problem initially had some constructive implication, because 
in [6] was suggested a special new kind of stochastic quantum networks to resolve 
the problem of universality. In the stochastic programmable quantum gate arrays [6] 
the size of a program is also only in two times bigger than the data, but the result 
of calculation is non-deterministic. The stochastic, probabilistic design described in 
many works [6, 9, 10, 11, 13] and will not be discussed here with detalls. 

The models with mcasurcmcnts and probabilitics makc actual using programmable 
networks with mixcd states and density matrices [12, 14], but it also deserves a separate 
account and not prcscnted here. 

In the present paper is considered the deterministic design with pure states and 
it is shown, that the problem with universality has rather formal meaning. From the 
one hand, only finite number of diffcrcnt operators are necessary for the universality 
in the approximate sense often used by default already in the earliest works about 
the quantum networks [1, 4, 15] and so it is possible to construct an universal (in 
the approximate sense) programmable quantum network with finite program register 
[16, 17, 18, 19]. 

From the other hand, even if to use the notion of the exact universality [20] with the 
infinite number of programs, it is possible to use the quantum computations with contin- 
uous variables [21] for generalization of considered models for the infinite-dimensional 
program register [18]. 

The models of programmable quantum networks with pure states considered in this 
paper have the direct analogue with the conditional quantum dynamics [22]. 
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In Sec. 2 is represented the formal theory of programmable quantum networks with pure 
states. Definitions are revisited in Sec. 2.1. Limitations due to laws of quantum me- 
chanics and thc conception of the deterministic and stochastic programmable quantum 
gate arrays are discussed in Sec. 2.2. Quantum networks presented in this paper are 
deterministic in such classiíication, but some specific points relevant to the stochastic 
design may be found in Sec. 2.3. The deterministic design is originated from the idea 
of conditional quantum dynamics recalled in Sec. 2.4. In present paper arc also mainly 
are used quantum systems with finite-dimensional Hilbert spaces, but a hybrid model 
with continuous quantum variables is also not omitted and may be found in Sec. 2.5. 
The design of a programmable Control-Shift network is represented íinally in Sec. 2.6 
and is used as a bàsic model for the rest of the paper. 

In the Sec. 3 are discussed more concrete models. The theory of universality in the 
quantum computations and control is recollected in Sec. 3.1. All nccessary Hamilto- 
nians for the programmable Control-Shift network with the universal set of quantum 
gates are constructed in Sec. 3.2. The relation with the theory of quantum cellular 
automata is demonstrated in Sec. 3.3. 

2 Formal theory of programmable quantum net- 
works 

2.1 Definit ion 

Let us return to the main equation describing a programmable quantum network 
Eq. (1). It is also may be written in more general form, when the state of the program 
register may be changed after the operation: 

\J: {\D)(g)\n)) ^ {un\D)) (g)\n'). {!') 

The scheme of a network for such transformation is depicted on Fig. 3. In the 
Eq. (1') is used expanded notation with the tensor product \D) ® |i7) often omitted 
for simplicity in expressions like \D)\n) or \D,n). 

Sometime it is convenient to exchange program and data registers and to write 

U: (|i7) ® \D)) ^ \n') ® {un\D)). (1") 

2.2 Limitations due to quantum laws 

The Eq. (1) and Eq. (1') was analyzed in [6], there such kind of networks was called the 
programmable quantum gate arrays. Sometimes they are also called quantum processors 
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Figure 3: Scheme of a programmable quantum network 

[8, 11, 12, 13, 14, 17, 18, 19]. In Eq. (1') the state of the program register |i7') after the 
operation should not depend on the state of the data register, because for two different 
data States it could be written 



U(|í?i)S 


D |i7)) = 




■o\n[) 


\J{\D2)^ 


D |i7)) = 


{un\D2)) S 





(3) 



The unitary operator U does not change the scalar product and it is possible to write 

{D,\D2) {n\n) = {D,\D2){n[\n'2). (4) 
1 

If States \Di) and \D2) are not orthogonal {Di\D2) 7^ 0, Eq. (4) is satisfied only for 
(Í7í|i7^) = 1 and so \n[) = |i7^). 

For (D1ID2) = such an argument does not work, because Eq. (4) is always true. 
It is really possible to resolve Eq. (3) for orthogonal data states, but in such a case for 
different states \n[) and lilg) due to linearity it may be written: 

U((a|Di) + (3\D2)) ® |i7)) = a{un\Di)) ® \n[) + í3{un\D2)) ® |i7^), 

so the state of program and data registers may be entangled and the equation does not 
have required form. Only for \n[) = lU!^) = |i7') the equation may be reduced to the 
proper form Eq. (1') for the arbitrary superposition of data states 

U + ^\D2)) ® |i7)) = {un{a\Di) + (3\D2))) ® |i7'). 

In the [6] also was found another important consequence of the given structure of 
Eq. (1') — the orthogonality of different program states. Really, let as consider two 
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diíferent programs |i7) and \S). It is possible to write 



U(|L>)(8)|i7)) = {un\D))^\n') 
\J{\D)®\S)) = {us\D))®\S'). 



(5) 




(6) 



In [6] was noted, that because only the term {D\unUs\D) in Eq. (6) depends on the 
state \D), it may be resolved only for {n\S) = = {n'\S') i.e., if all different programs 
correspond to orthogonal states. So the orthogonality condition Eq. (2) mentioned in 
Sec. 1 is provcd. 

• Note: It should bc mentioned, that the question about behavior of tcrms hke {D\u\D) 
maybe denotes additional discussions, but it is above a scope of presented paper. Say, it is 
possible to consider a real-valued analogue of the qubit [23, 24], a "reòíí." It is a vector \R) in 
the two-dimensional real vector space (plane). The operators u now correspond to rotations 
of the plane. In such a case {R\u\R) = cos(<^) does not depend on the state of the rebit (here 
is the angle of rotation). 

Due to the orthogonality property Eq. (2) for a program register in Eq. (1), the 
dimension of the Hilbert space is equivalent to the number of different operators Un- 
For the exact universality we formally must have possibility to apply the iníinite number 
of different unitary operators u. It was considered in [6] as some disadvantage and it 
was suggested idea of "probabilistic" (non-deterministic) programmable quantum gate 
arrays. 

2.3 Non-unitary linear operators and non-deterministic net- 



It should be mentioned, that the no-go result about the exactly universal deterministic 
programmable quantum network may be considered as a more subtle limitation in 
comparison with a famous quantum no-cloning theorem [7]. The no-cloning theorem 
may be derived from the linearity of the quantum mechanics, but the consider ation 
above also uses the unitarity. It is essential, because exist linear, but non-unitary 
operators satisfying Eq. (1') with the finite-dimensional program register. 

If a data register is described by some state vector \D) & 7í with components 
Di, then the minimal program register must contain the matrix of coefíicients for the 
operator A^j and so in the simplest case such a register in two times bigger than the 
data one \A) eH^H. Here due to the usual law of composition of quantum systems 
for A^-dimensional data register such a program register is A'^^-dimensional. 

Let \A) = Aij\i)\j) is the state of the program encoding the operator A = 
'^Aij\i){j\, viz \ A) = |j)A|j) and M is the linear non-unitary operator defined on 
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the basis as 

M: \k)\t)\j)^ó,,\t)\Q)\Q). (7) 
Then a specific kind of Eq. (1') is vàlid 

M(|^)|A)) = |A^)|0,0), (8) 

because {Aip)i = J2j ^iji^j = Xljfc ^ij^jk'ipk- In iact, the M is the operator of the matrix 
multipUcation rewritten in a specific way. 

The idea has some relation with the non-deterministic design [6] . Definition of the 
state \U) representing an unitary operator U may be simply changed, to ensure the 
unit norm 

i 

For the qubit (two-dimensional Hilbert space) it coincides with the definition in [6]. 
It is possible to rewrite \ip)\U) using the Bell basis $± = ^(|00) ± |11)), *± = 
^(|01) ± |10)) for the first two qubits [6] 

1^)1^) = + \^^)UaM + + |$-)L/o-,|^)). (10) 

With minimal modification such expression may be converted to the expansion of 
some unitary operator G = L + R, there L is the non-unitary linear "programming" 
operator and R is the residue R = G — L 

G{mu))-{umM+R{mu)). (n) 

So the non-unitary linear operator at least formally presents in the model of the non- 
deterministic programmable quantum gate array, but further discussion on the question 
is outside of the scope of this paper. 

2.4 Conditional quantum dynamics 

Let us return to the programmable quantum networks denoted in [6] as deterministic. 
The result about the orthogonahty of different program states also has some positive 
implications and simplifies description. It is possible to choose the states corresponding 
to different programs as a new basis and to write instead of {11^) simply \k). Structure 
of the operator U Eq. (1) becomes quite clear in such notation, it is the conditional 
quantum dynamics [15] described even earlier, than the programmable quantum gate 
array s mentioned above. 
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Such an unitary operator U for orthogonal \k) from m-dimensional Hilbert space 
may be written [15] 



m—l 




U = ^ \k){k\®Uk, \JÍ\k)\D)) = \k){uk\D)). (12) 

k=0 

In Eq. (12) program and data registers are exchanged in comparison with definition 
Eq. (1) or Eq. (1') and it corresponds to the alternative notation Eq. (1"). Such an 
order is used further in the paper for convenience (see Note on the page 11). 

A simplest example — is the controlled-NOT, c-NOT gate. 
It was used already in earliest papers about quantum com- 
puters [25] and also called the measurement gate [1]. It is Cl Cl 
the quantum version of a classical reversible gate with two 
input and two output states. The c-NOT applies NOT to a 

second state if and only if the first state is 1 (TRUE), i.e., for — G^Ò 

Boolean variables a, b may be described as (a, b) i— > (a, b@a), 
where "@" is addition modulo 2, 0@0 = 1@1 = 0, 
0@1 = 1@0 = 1. 

The quantum mechanical version is straightforward 

c-NOT,^: |a)|6) ^ \a)\b@a). (13) 

Two possible instances of the c-NOT gate are depicted on Fig. 4. They also may be 
considcred as representations of two different orders of program and data registers used 
in Eq. (12) and Eq. (1). 




Figure 4: Schcmcs of controlled-NOT quantum gates, a) The first qubit — is control 
\a)\b) 1-^ |a)|ò0 a). b) The second qubit — is control \a)\b) ^ \b@ a)\b). 

It is useful to consider the simple case to show principle of matrix representation 
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of such operators. The action of the c-NOT gate for the basis may be written as 



|0)|0) 

|o)|i) 

|1)|0) 

11)11) 



|0)|0) 
|0)|1) 

|1)|0) 



so the c-NOT gate simply exchanges two last vectors of the basis |10) 
matrix of the operation may be written as 



(14) 



and the 



c-NOT,, 



/ 1 








1 



|00) 
|01) 
|10) 



(15) 



UnUke the classical case, the difïerence between c-NOT^s and c-NOTa, quantum 
gates Fig. 4 is rather formal — it is enough to change bases for the both systems to 
convert one gate to another and so both qubits formaUy are equal in such operation 
[27, 28]. The new basis is 



l+) = ^(|0) + |l)), 



11))- 



(16) 



In the new |±) basis c-NOT^a Eq. (14) is rewritten 



= è (|oo)+|oi)+|io)+|ii)) (|oo)+|oi)+|io)+|ii)) 
= è (|oo)-|oi)+|io)-|ii)) (|oo)-|oi)-|io)+|ii)) 

= i (|00) + |01)-|10)-|11)) ^1 (|00) + |01)-|10)-|11)) 
|-)|-) = i (|00)-|01)-|10) + |11)) ^1 (|00)-|01) + |10)-|11)) 

The Eq. (17) show, that a first element is swapped if and only if the second one is 
|— ), it is just C-NOT21 written in the new bases |±)|±). So there is no clear distinction 
between the control and the controUed system. It is shown further, that such situation 
may be not vàlid for more difficult cases. 

Similarly with controlled-NOT it is possible for an arbitrary gate C/ = ( ) to 
define the controlled-Z7 gate sometimes denoted as Ai(Z7) [26] 



-)|-) 



(17) 



+ 



u 



/ 1 \ 

10 

Uoo Uoi 

V C/10 Uu ) 



|00) 
|01) 
|10) 

lll) 



(18) 
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Using the same principle, the sum used in Eq. (12) may be written as the block- 
diagonal matrix 



U = 



V 



Uq 











\ 



Ur 



(19) 



-1 / 



• Note: For the m-dimensional program register and the n-dimensional data register such U 
is mn X mn matrix with only nonzero elements are n x n blocks with matrices on diagonal 
of U. The understanding notation Eq. (19) may be used if the program register is the first 
system like in Eq. (1") and Eq. (12). 

Due to the orthogonality condition Eq. (2) the dimension of the program register 
often much bigger noi maybe even infinite and so, unlike the example with c-NOT 
gate, two systems used for data and program registers in general are not equivalent. 



2.5 Infinite-dimensional program register 

It was already mentioncd, that for the universality in the exact sense the dimension 
of a program register should be infinite [6]. Let us consider such a case. Such a 
programmable network uses both continuous and discrete quantum variables for the 
program and the data registers respectively [18] and so may be considered as an example 
of the hybrid quantum network [29]. 

The system is hybrid also in other meaning [30], i.e., it is the possible approach to 
the unified description of a quantum system controUed by some analogue parameters. It 
make some bridge with usual (non-programmable) quantum networks based on pseudo- 
classical description of "tuneable" gates depicted schematically on Fig. Ib. 

For transition to the continuous variables it is possible to change the sum in Eq. (12) 
to the integral and write formally [18] 

^ = ji\q){q\®n{q))dq, [J{\q)\D)) = \q){uiq)\D)). (20) 

The Eq. (20) describes the family of gates (matrices) u{q) parameterized by some 
continuous variable q. Elements of the basis in infinite-dimensional Hilbert space are 

denoted as \q). 

A Standard example — is the space of functions ip{x) on a line and Dirac delta 
functions as the basis 

\q) = S{x-q), {q\i^{x))=i:{q). (21) 

The sign ® of tensor product may be skipped in the expression like Eq. (20), because 
for such a hybrid case with continuous and discrete quantum variables there is simple 
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interpretation. The tensor product of the space of functions on a line and some finite- 
dimensional vector space may be represented as a space of multi-component functions 
on the hne with vàlues in the vector space. Say the both program and data registers 
may be encoded in a wave vector of one particle. 

The space of hnear operators in such a case may be represented as a space of 
matrix-valued operators on the hne. For example Eq. (20) may be expanded as 



\j{'ilj{x)\D)) = / 6{x - q)'^{q)u{q)\D)dq = 'ilj{x)u{x)\D). 



(22) 



Using the example with one particle, Eq. (22) may be considering as a process with 
"twisting along x" of the wave function |^E'(a;)) depicted on Fig. 5. 





|*(x)) = ipix)\D) 



U|^(x)) = ip{x) u{x)\D) 



Figure 5: A scheme of the process Eq. (22) with a "distributed" qubit. 
It is also possible to use the momentum basis of the periòdic functions 

\p) = exp{2nipx). (23) 
With the new basis it is possible to write 

= 1 {\p}(í^\u{p))dp, \J{\p}\D)) = \7){uÍp)\D)), (24) 

there is also taken into account the possibility of change of the momentum p —>■ p'{p) 
(the program register) similarly with Eq. (1'). 

More detalls may be found in [18], let us only consider a visual example with "a 
scattering process" Fig. 6, there are used two different systems. 

Before the scattering (Fig. 6, time ti) the first system was described by a 
"monochrome" wave function exp(ifca;) and the second one — is a localized qubit 
\D). After the "inelastic" scattering process (Fig. 6, time Í2) the first system may have 
a state with other k' = k'{k) and the qubit changed a state to u{k)\D). 

Such a scattering process maybe too formal and abstract, but it provides the under- 
standing analogue between the programmable quantum network and a more traditional 
design with a quantum system controlled by some external devices like làsers. 
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Figure 6: lUustration of Eq. (24) with scattering process. 

It should be mentioned, that the considered model with continuous variables also 
provides some link with the stochastic programmable quantum gate arrays. It was 
already mentioned, that such a model produces a correct answer with some probability. 
For A^-dimensional data register the probability is i.e., 2^" for a system with n 

qubits [6]. For some particular case of encoding U{1) operations, it is possible to 
make the probability higher for the bigger program register with approach to the unit 
probability of success in the infinite limit [9, 10]. 

On the other hand, for the deterministic design it is also possible to consider the 
infinite-dimensional program register as a limit of finite-dimensional one. Sequences 
of finite-dimensional networks used for such limits have quite different properties for 
stochastic and deterministic approach. For the first one we may apply an arbitrary 
operation U{1), similar with rotations on arbitrary angles 2n(j), but it succeeds only 
with some probability increasing with approach to the infinite limit. 

For the second one, the operation is always successful, but we may apply only finite 
number of different operations U{1), i.e., rotations on some fixed angles 27rk/N, and 
for the infinite limit the rotations cover the full circle. Despite of such difference for the 
finite case, the infinite limits for both deterministic and stochastic networks for U{1) 
operation are essentially the same. A rather technical proof of the interesting fact may 
be found in [18] and is not reproduced here. 

2.6 Programmable quantum networks with three buses 

It should be mentioned, that Eq. (1) and all derivative equations considered here are 
still do not have the complete analogue with the idea of usual (von Neumann) computer 
architecture, because formally they are describe only one step of a program. 

From the one hand, such a picture sometimes may be appropriate for the theory of 
quantum computing. On the other one, it is very limiting, because even in the Standard 
definition of the universal set of quantum gates, it is used the arbitrary composition of 
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the gates from the set [1, 4, 15, 20, 26]. 

Such an idea of the universality may be implemented also for the programmable 
quantum networks with pure states [16, 17, 18, 19] discussed in this paper. The scheme 
Fig. 7 uses two operators and three buses. 
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Figure 7: The programmable quantum network (processor) with three buses 

The C, Control — is the aheady considered operator Eq. (1") of one computational 
step with second (control) and third (quantum data) buses. The S, Shift operator must 
change the state \k) of a control bus after each step using first (program) and second 
buses. 

The buses also may be called pseudo-classical, intermediate and quantum, because 
the quantum bus may contain arbitrary superposition of states, the intermediate bus 
should use only orthogonal states, but it is still may not be considered by an entirely 
classical way, because it is linked via C operator with the quantum bus and the specific 
character of such a design was already discussed above. But already for description of 
the pseudo-classical bus and the operator S it is at least formally enough to use the 
theory of reversible classical gates. 

In principle, S may correspond to an arbitrary reversible classical program, but here 
is convenient to consider the simplest case with the cyclic shift. Let us suggest, that a 
data register is described by m-qubits, there is a finite set with N = 2'^ quantum gates, 
and for our purposes it is necessary to apply a sequence with L gates from the given 
set. For example we may consider the finite set of universal gates and the task of the 
approximation of an arbitrary gate with the necessary precision using up to L gates. 

For the design with the cyclic shift operator S a program register must have size 
{L — l)n qubits and acts on Ln qubits of program and control registers as 

S: \kL,kL-i,...,k2)\ki) ^ fc^, . . . , /cg) j/cs). (25) 
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Action of C was already defined in Eq. (12) 

C: \k)\D)^\k){uk\D)). (120 
and so the composition of C and S may be written 

Sc(\kL,kL-i,...,k2)\h)\D)^ = \ki,kL,...,k3)\k2){uk\D)). (26) 

It is only one step. Let us denote the state of program and control registers as \K) = 
|/cl, ^2)1^1)- For L steps \K) returns to the initial state and so it is possible 

to write [16, 17, 18, 19] 



(SC)^: \K)\D)^\K){uL---U2Ui\D)) 



(27) 



Here the program bus should be rather compared with the cyclic read-only memory 
(ROM). A scheme of the programmable "Control-Shift" network with the cyclic shift 
register is depicted on Fig. 8. 



Program (cy^lfó ^ÜQM) 




Figure 8: A scheme of Control-Shift network 
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3 Models of programmable quantum networks 



3.1 Universal sets of quantum gates 

In the theory represented above were used rather abstract models. Standard relation 
between the scheme of a quantum gate and the evolution of the quantum system is 
depicted on Fig. 9. The input and output of the network — is really the same system 
at two different moments of time and the gate represents change of the state of the 
system during given period due to some interactions. 




Figure 9: A scheme of a gate and evolution of a quantum system 

Idea of decomposition using the universal set of quantum gates was already few 
times mentioned in the present paper, e.g., it is relevant to the model of the Control- 
Shift network described by Eq. (27). Let us recali the bàsic principies of the theory of 
the universality in quantum computations and control. 

There is so-called infinitesimal approach with the Lie algebras [15, 31, 32]. It 
is convenient also due to the direct relation with Hamiltonians of quantum systems. 
For the system with the constant Hamiltonian H the evolution during the time t is 
described by the unitary operator (gate) 

U = U{t) =exp{-iHT). (28) 

If two Hamiltonians Hi, H2 correspond to gates Ui, U2, then 

t7it/2 = exp(-ïiíir)exp(-ïiÍ2r) =exp(-z(iíi + iÍ2)r) +0(r2), (29) 
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there O(r^) is an error of order and so for small r the sum of Hamiltonians corre- 
sponds to the composition of the gates. It is also quite Standard to use an expression 
for commutators [15, 31, 32] 

U.U^U^'U^^ = exp(-(iïiiÏ2 - H^H^y) + 0{t^) (30) 

and so such a product with four terms is approximately equal to the action of a gate 
with a Hamiltonian — i[13'i, ^"2] and a parameter r^, ie., if to consider t — y/t the 
Eq. (30) has prccision 0{t^-^). 

Due to Eq. (29) and Eq. (30), which vahd for infinitesimal vàlues r ^ 0, it is possible 
to formulatc a condition of universality using Hamiltonians [15, 31, 32]. The set of 
Hamiltonians Hk corresponds to the universal set of gates, if it is possible to produce 
an arbitrary Hamiltonian using linear combinations of Hf. and they commutators of 
any order, i.e., i[Hj,i[Hk, . . .]\. 

Such dcfinition is convenient for a test of universality, but for concrete tasks it may 
produce a problcm due to equations for commutators. Already for Eq. (30) with 
it is not clear, how to choose a parameter r for the generation of a gate with a good 
precision, and the problem even worst for commutators of k-th order with t''. 

There is a way to produce an expression with the íirst degree of r instcad of 
in Eq. (30) and without error at all if to use special choice of Hamiltonians [33]. Let 
us consider a systcm with n-qubits. As a basis in the space of Hermitian matrices 
(Hamiltonians) may be used 4" different tensor products of four Pauli matrices 

Hj^(Tj,®(Tj^0---0(Tj^, jfe = 0, ...,3 (31) 

where j = (ji, j2, • • • ,jn) is notation for multi-index and 




It is also widely used the alternative notation: 1, cr^, (Ty, cr^ respectively. 
All Hamiltonians Eq. (31) have properties 

iïji/k = iiïkiïj, ííj = 1. (33) 

Unity of the square in Eq. (33) ensures a simple expression for the exponent 

exp(i0i3"k) = cos(0)l + ism{(j))Hk (34) 

and using the identity 

exp{A) exp{B) exp(— A) = exp(exp(A)Bexp(— A)) 
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it is possible to find the precise expression for the exponent of commutator 



Unlike Eq. (30) the Eq. (35) contains the first degree of r and for such a method of 
generation of arbitrary gates there are only errors related with Eq. (29), which have 
order Oii?). 

It is well-known that one- and two-gates are enough for universahty [15, 26, 31]. 
For the work with Hamiltonians hke Eq. (31) with properties Eq. (33) it is convenient 
to use an universal set with 2n + 1 elements [33] 

Zq = (7^^^ = CTs (g) q-Q íg) ■ ■ ■ (g) Q-q ^ , (36a) 

n-l 

(2) 

Zi = Cr^ ^ = (To (g) (73 (g) CTo (g) • • ■ (g) CTo , (36b) 



n-2 

y ^ 



= (Ti = (To (g) ■ • ■ (g) CTq (g) (Ti (g) (To ® ■ ■ ■ (g) (To, , (36c) 

Dk — crf^^^a^^^'^'^ = (To <g) • • • <g) <To (g) íts (g) (T3 (g) (To (g) • • • (g) (Tq . (36d) 



n-fe-2 



Really the set of gates diífers from suggested in [33] on nonessential change of the 
basis ((Ti (T3). After such exchange all two-qubit gates have diagonal form due to 
Eq. (36d). It should be reminded, that the gates are generated using Eq. (28) with 
given Hamiltonians Eq. (36). It may prevent some confusion, because all operators in 
form Eq. (31) are not only Hermitian, but also are unitary and so they are widely used 
as gates in many works. 

There is also alternative approach to the universahty, then instead of consideration 
of Hamiltonians and infinitesimal (small) transformations, it is considered a question, 
how to decompose precisely the given unitary matrix (with 4" parameters for n qubits) 
on product of one and two-qubit gates. Standard choice is: all one-qubit gates to- 
gether with c-NOT gatc [34, 35, 36], e.g., in [36] is suggested an algorithm for the 
decomposition with 4" one-qubit gates and 4" — 2" c-NOT gates. 

It should be mentioned, that c-N0Tl2 gate also may be transformed to diagonal 
form, if to change the basis of the second qubit to |±) Eq. (16). Such transformation 
is described by the Hadamard matrix H = ^ { \ li) and already was used above on 
page 10, there it was applied to both qubits for the transition between C-NOT12 and 
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C-NOT21 gates. After such diagonalization the Hamiltonian Eq. (15) may be written as 



d-N0T = Ai(í^3)- 



/ 1 









-1 



(37) 



Really, any control Ied-Í7 gate Eq. (18) may be transformed to diagonal form by some 
change of a basis, ie., using one-qubit operations with the second qubit. It is simply 
a basis there U itself is diagonal. 



3.2 Hamiltonians of element ary programmable gates 
3.2.1 Control gates 

Let us consider models of bàsic quantum gates, necessary for implementation of the 
universal programmable networks. In Sec. 3.1 it was shown, that it is enough to use 
"powers" of Pauh matrices 



exp(icrfer) = cos(t) + isin(T)<Tfe 



(38) 



together with some simple diagonal two-qubit gate. In the programmable quantum 
networks are used controUed gates and so bàsic elements are two- and three-gates. 

The Hamiltonian of controUed versions of a gate may be found using the simple 
expression 



1 








(39) 



where 0, 1 are zero and identity matrices respectively. So Hamiltonians of two- gates 
for the controUed one-qubit operations used in the universal set Eq. (36) may be rep- 
resented as 
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(40) 



The three-qubits Hamiltonian for the controUed diagonal gates Eq. (36d) is 
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Let us show, how to use instead of Eq. (41) only two-qubit Hamiltonians. 

= l(.(^).f-c.i^)c.fc.f) = 1(^(^)^(3) (42) 

If to consider thc three-qubit gate generated by h^, due to the decomposition Eq. (42) 
and Eq. (35) it is possible to write 

U^t) = eM^h^or) = exp(^^^f ) exp(-^^^Wa(2)^f ) 

= exp(z^crf ) exp(z^a(Vf ) exp(-z^^W^(3) ) exp(-z^^WcTf ). (43) 

Let us consider also a Hamiltonian with two control qubits for implementation of 
all three Pauli matrices 
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(44) 



It is the simple example of a Hamiltonian for the programmable network with four 
diíferent "programs" and one data qubit. If first two qubits have the state |0)|0), the 
data qubit is not changed, but if they have state |0)|1), |1)|0), or |1)|1) during a period 
T, then the data qubit is changed as \D) i— > exp{icrkT)\D) with k = 1,2,3 respectively. 

From the one hand, the example shows a specific kind of the exactly universal 
programmable quantum gate satisfying Eq. (12) with the finite-dimensional control 
register. Here the continuous parameter r is used to implement the infinite number 
of different programs. From the other hand, for a programmable quantum network 
it is more appropriate to fix some period Ar and use conseqüent application of gates 
Uk = exp(icrfcAr) via the three-buses design discussed in Sec. 2.6. 

Here is not suggested, that Ar should be infinitesimally small. Using the infinites- 
imal parameter is not the only way of approach to a continuous limit. Most Hamil- 
tonians used here generate gates with a simple periòdic behavior due to Eq. (34) and 
if to choose Ar as an irrational múltiple of tt, it is possible to approximate any real 
parameter with an arbitrary precision. Such an idea really was used already in the 
earliest works about the theory of universal quantum computations [1, 4, 15]. 
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Using qubits for the control is convenient, but is not necessary. E.g., it is enough for 
universality to use only two Pauli matrices, and so instead of two qubits it is possible 
to use one qutrit for '3— >2' control, i.e., a quantum system with three states |0),|1) and 
|2) and to write a Hamiltonian like 
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(45) 



The very rough scheme of sucli a system is depicted on Fig. 10. Here a qubit — is 
a system with spin-1/2, and a qutrit is represented as a quantum system distributed 
in a potential of a triple quantum dot "molecule." It is suggested, that in the state 
|0) the qutrit is not interact with the qubit, but the states |1) and |2) already affect 
on the spin system and produce some evolution of the qubit like \D) i— >• exp{i(TkT)\D), 
k = l,2. 




.4.. '^^^^ ^^^^^^ 

Figure 10: A scheme of six bàsic states for '3— >2' control Eq. (45). The qubit is depicted 
as a system with spin and the qutrit as a "triple quantum dot." 

Let us return to the controlled two-qubit gates Eq. (41). It is convenient for further 
applications to place the control qubit between the two controlled ones and then, 
rewriting Eq. (41), it is possible to consider a diagonal Hamiltonian with three energy 

E = : |0)|0)|0),|0)|0)|1),|1)|0)|0),|1)|0)|1) 

E = -AE : |1)|1)|0),|0)|1)|1) . (46) 

E = +AE : |0)|1)|0),|1)|1)|1) 

The scheme has visual interpretation: if the control qubit is in the state |0), the two 
data qubits are not interact, but if the control qubit has the state |1), then the energy 
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of the system is bigger for the data qubits in same states and smaller for different 
ones. So the program qubit may be considered as a "control switch" of the interaction 
between the two data qubits. An illustrative, but a rather naïve scheme with a double 
quantum dot for the control of two spin systems is depicted on Fig. 11. 




Figure 11: A scheme of energy levels for the diagonal Hamiltonian of control Eq. (46) 

Some discussion on physical systems appropriate for such a "switching" purposes 
may be found also in papers related with models of the "global" quantum computing 
(with "always-on interactions" ) [37, 38, 39, 40], because such approach is very close to 
the idea of programmable quantum networks. 

Finally, a programmable network with such a kind of gates for one- and two-qubit 
operations is depicted on Fig. 12. Here a data register is presented as an array of spin- 
half systems and a control register is consisting of double and triple quantum dots. 
Double quantum dots used for two-qubits operation here are intermittent with data 
qubits, and qutrits are situated above. 

The controlled one-qubit quantum gate for such a register may be described by 
Hamiltonians Eq. (40) and Eq. (45) 

K= |0)(0|®l + |l)(l|®cr, + |2)(2|®cr,. (45') 

It is enough simply to include number k of controlled qubit for such Hamiltonian h'^^''\ 
On the considered scheme Fig. 12 for data qubits the number is A; = 2j ' + 1. 
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Figure 12: A visual scheme of a programmable quantum register 

Hamiltonians of controlled two-gates devote a bit more detailed consideration, be- 
cause one data qubit is controlled by two different gates. Let us rewrite Eq. (42) for 
the control qubits between the two data qubits 

hl = a.®r,^,)®CT. = l {a^pcrf - ^IVf ) . (47) 

The analogue of decomposition Eq. (43) is 

Ul{r) = eM^h^r) = exp(.^a«^f ) exp(-^Ia«^f ) 

= exp(.^^(Vf)) exp(z^^«ai^)) exp(-.I^(^)^f ) exp(-^^^(V(2)). (48) 

For an array with qubits Hamiltonians Eq. (47) may be written as 

The Hamiltonians Eq. (49) commute for different k and so despite of the overlap, again 
there is no an essential difference between the control of array Fig. 12 and the initial 
model with only three systems. 

3.2.2 Shift gates 

The models above represent only C (control) part of the Control-Shift network discussed 
earlier in Sec. 2.6. It is necessary to supply new and new indexes \k) for the control 
register and if control and data buses may be considered as some arrays along a given 
axis (x) Fig. 12, it is possible to arrange a program bus as some bar along an orthogonal 
axis (y) with a separate line (ID array) for each control gate and to produce the two- 
dimensional array Fig. 13. 

In such 2D structure, each control element is receiving indexes from only one line, 
i.e., ID array along orthogonal axis {y, program). Let us consider necessary operations 
with the array. 
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Figure 13: A two-dimensional programmable quantum network 



The simplest idea — is to implement the cyclic shift used Eq. (25) for each such 
ï/-array 

S: \kL)\kL-l) ■ ■ ■ \h)\h) ^ \h)\kL)\kL-l) ■ ■ ■ \k2). (25') 

Such a method produces vahd reahzation of a programmable Control-Shift network, but 
there is a difíicuhy with implement ation of the Hamiltonian for Eq. (25'), because it is 
not a /oca/ operation, z.e., it acts on all L systems. 

On the other hand all Hamiltonians for the control above were local with three 
systems or less. It is also possible to write local anàlogues of the Shift operation. The 
simple way — is to consider two different operations: for one step are exchanged all 
pairs |A;2i+i) and 1^2^+2) and for next step — pairs |A;2j+2) and \k2j+z)- 

Si : ■ ■ ■ |/C2j+l)|/C2i+2) ■ ■ ■ ■ ■ ■ \k2j+2)\k2j+l) ■ ■ ■ /ggN 

S2 : ■ ■ ■ I k2j+2) I k2j+z) ■ ■ ■ ■ ■ ■ I fcaj+s) | ^2^+2) ■ • • 

So each step is performed using only the exchange (SWAP) operation with two 
systems. For qubits, qutrits, and quantum systems with an arbitrary number of states 
( "qurfits" ) the SWAP gate is defined on basis states as 

SWAP: |A;)|0 ^^ lOI'^)- (51) 
For qubits it may be also implemented using three c-NOT gates [25] 



SWAP 



/ 1 \ 

10 

10 

\ 1 / 



c-NOTi^c-NOT^iC-NOT,^. (52) 
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If to use the alternating sequence of such operations, states with odd numbers 
like \ki) are permanently shifted in one direction and states with even numbers — in 
opposite one. So it is possible to encode a necessary sequence if to use only the odd 
states and to set other to zero 



Here is suggested as usual, that the zero index corresponds to the identity operator. 
3.3 Quantum cellular automata 

The present paper is not devoted exclusively to the theory and applications of quantum 
cellular automata (QCA), but it is reasonable to discuss briefly the tòpic, because QCA 
models are widely used in many works devoted to related problems. It was already 
briefly mentioned the global quantum computing [37, 38, 39, 40]. Recent works on 
reversible quantum cellular automata [41, 42, 43] make possible to talk about even 
more direct relation between the model of universal programmable quantum networks 
with pure states and quantum cellular automata. 

The just considered model Fig. 13 illustrates, that the Control-Shift design has very 
close relation with cellular automata, but there is specific subtleties for transition to 
the quantum case [41, 42]. 

Let us consider for example a spin lattice. It could be used formally for represen- 
tation of a cellular automaton with two vàlues {0, 1} encoded by two bàsic states | j). 



(53) 



i) Fig. 14. 




Figure 14: A lattice with qubits representing a cellular automaton 
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On the other hand, a general state of such a lattice is not necessary may be rep- 
resented as a product state Fig. 15, then it may be simply treated as a collection of 
celis in different states. The general entangled state is defined as a sum Fig. 16 with 
complex coefScients on different basis configurations Fig. 14. 



Figure 16: A sum on different configurations representing general entangled state 

This specific property of a quantum system does not permit to talk about a state of 
a single celi, or a block of celis as about a vector in the Hilbert space, and it produces 
specific difficulties for local description of cellular automata. 

On the other hand, it was already mentioned, that the program bus may be formally 
described using reversible classical computations and the theory of reversible cellular 
automata has a long history and is appropriate for such a purpose [44, 45, 46]. 




Figure 15: A lattice with qubits described by the product state 
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One specific method used in the theory of reversible cellular automata [45] is the 
Margolus partition with two different operations during different time steps. The sim- 
plest example for ID automata just corresponds to Eq. (50) used for local implemen- 
tation of the shift operation. 

The alternating sequence of S and C operations in Eq. (27) used above in pro- 
grammable Control-Shift networks in Sec. 2.6 also becomes a quantum analogue of the 
Margolus partition, if to consider 2D design of quantum gates described below Fig. 13 
3jS cl quantum cellular automaton. 

The design considered here maybe still does not look like "true" QCA, because it 
is not regular enough, e.g., there are 2D array for program with two different kinds of 
systems and ID array of data. On the other hand, cellular automaton with different 
kinds of celis is equivalent to cellular automaton with only one kind with additional 
fiag. 

It is also useful to consider 2D array of data [43]. We just need to spread the whole 
programmable register Fig. 12 along program bus Fig. 17 (instead of spreading only 
control elements of the register depicted earlier on Fig. 13). 




Figure 17: A quantum network with two-dimensional data and program registers 

Such array may be considered as many independent copies of the programmable 
register. To produce initial design Fig. 13 from such a regular QCA, it is possible 
formally to introduce an additional fiag-qubit and to modify the control register in 
such a way, that it changes state of data only if the fiag is |1). 

On the other hand, it is possible to save the initial structure and to let all program 
registers evolve in paral·lel. Implementation of the cyclic shift may be represented by 
cylindrical QCA Fig. 18. 

Here we have three bàsic elements: the qubits of quantum data register (s), triple 
quantum dots for implementation of one-qubit gates (below a qubit on Fig. 12, Fig. 17, 
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Figure 18: A programmable quantum network as cylindrical QCA 



etc), and double quantum dots for control of two-qubit gates (between each two qubits 
on Fig. 12, it should be mentioned, that on Fig. 17 the registers are arranged along the 
axis of the cylinder). 

The QCA programmable network uses two Standard steps C and S of Control-Shift 
network, but they are subdivided on smaller steps. The analysis of operation S is 
simpler due to direct relation with reversible classical computations already mentioned 
above. 

The S(hift) steps move data around cylinder and may be implemented using two 
sub-steps Si,S2 Eq. (50). It was already mentioned, that the indexes in such a case 
should be intermittent with zeros Eq. (53). 

Each control line may be considered independently and so S step is described by 
the expression Eq. (50). For implementation of the universal set of gates like Eq. (36) 
two different kinds of lines are necessary: for the control of two-qubit operations are 
used only two indexes {0, 1} and for controlled one-qubit gates (at least for two first 
qubits) are used three indexes {0, 1, 2}. 

The lines with indexes for one-qubit gates are alternating with lines for two-qubit 
gates. It is important to mention, that all two-qubit gates from Eq. (36d) are commut- 
ing and so may be applied on the same step. The one-qubit gates on different lines are 
also commuting. On the other hand, gates Eq. (36c) are not commuting with gates 
Eq. (36d) for adjoint lines, and so it is necessary to apply such gates on different steps. 
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It requires a special arrangement of nonzero indexes on adjoint Unes. 

Let us suggest that such arrangement condition is satisfied to prevent consideration 
of non-commuting Hamiltonians and consider C(ontrol) step. 

One-qubit operations may be defined locally for each pair with the control qutrit 
and the data qubit. The Hamiltonian of control on such a pair was already discussed 
above Eq. (45') on the page 22. The Hamiltonians of controllcd two-qubit gates again 
coincide with expressions for single controUed register Fig. 12 and also were considered 
earlier Eq. (49). 

It was already mentioned, that the Hamiltonians Eq. (49) for control commute for 
different sites. The commutativity is important for quantum cellular automata models 
[41, 42]. Let us suggest, that an entire Hamiltonian of QCA may be exprcssed as a 
sum of all local Hamiltonians corresponding to local celis. Only for commutative oper- 
ators such a method produces appropriate relation between global and local transition 
functions 

exp(z J2 h'^^r) = eMih'''M, C^giobai = H ^^1-1 (54) 

k,j k,j k,j 

The Margolus partition used in the theory of reversible QCA [41] just uses two sets 
of commuting Hamiltonians. In the model considered here may be found even four 
such sets. The Margolus partition used for S(hift) operation is corresponding to the 
two operators Si, S2 Eq. (50) and already was discussed above, but C(ontrol) operation 
also contains two sets of Hamiltonians corresponding to one- and two-gates 

kJ kJ 

Let us introduce two operators 

Cl = exp(iií^Ar), C2 = expiiH^Ar), (56) 

with Hamiltonians of one- and two-gates respectively and Ar is the fixed irrational 
múltiple of 71, then the programmable quantum network based on such QCA may be 
described by the periòdic sequence of operators like 

Uvi = C1C2S1 C1C2S2, Uiv = CiSi C2S2, (57) 

there the second operator Uiv with only four terms is using alternating ('checker-board') 
arrangement, when after each step half of indexes are zeros and so only one operator 
between Ci, C2 produces nontrivial result. It should be mentioned, that both operators 
Eq. (57) correspond to two steps of a program and so for the cyclical design with the 
perimeter 2L, it is necessary repeat such operators L times. 
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4 Conclusion 



In this paper was considered the theory of programmable quantum networks. They may 
be considered as quite appropriate models of future quantum processors. Such devices 
would not require the macroscopic equipment traditional for modern experimental 
research in the area of quantum Information processing, because it is suggested to 
use quantum systems for the program, control and data. 

Such architecture is also interesting from the pure theoretical point of view. It let us 
use the unified description of all processes necessary for functionality of such quantum 
devices and reduccs amount of problems related with consideration of transition from 
classical to quantum domain. 

It should be mentioned, that some methods used in this paper also may be found 
in the more general theory about relations between quantum and classical pictures 
[27, 47] and so application for the quantum Information science is quite justified. 

In the present paper wcrc uscd programmable quantum networks with pure states 
(sometime also called 'deterministic') and it produccs additional clarification and sim- 
plification. It is also actual because 'probabilistic' programmable quantum gate arrays 
already are presented quite completely in other works. 

The theory of programmable quantum networks is still in the state of development 
and some new interesting branches may appear. For example universal programmable 
quantum computers by definition may perform any transformations with data, but 
it corresponds to arbitrary manipulations with a quantum systcm. Duc to such a 
principle emphasized already in earliest works [4] there is no big diffcrcncc between an 
universal programmable quantum network and a quantum robot [48] with wide range 
of possible opcrations. 

Say, it is quite reasonable to consider the famous Wigner's question about possibility 
of "sclf-rcproducing units" using the theory of universal programmable networks. In 
fact, some preliminary analysis was already performed using both probabilistic [50] 
and deterministic [51] approach, but the intriguing problem most likely devotes further 
research. 
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